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Abstract 

This paper addresses an interval analysis based study that is applied to the design and the comparison of 3-DOF 
parallel kinematic machines. Two design criteria are used, (i) a regular workspace shape and, (ii) a kinetostatic 
performance index that needs to be as homogeneous as possible throughout the workspace. The interval analysis 
based method takes these two criteria into account: on the basis of prescribed kinetostatic performances, the 
workspace is analysed to find out the largest regular dextrous workspace enclosed in the Cartesian workspace. 
An algorithm describing this method is introduced. Two 3-DOF translational parallel mechanisms designed for 
machining applications are compared using this method. The first machine features three fixed linear joints 
which arc mounted orthogonally and the second one features three linear joints which arc mounted in parallel. 
In both cases, the mobile platform moves in the Cartesian x — y — z space with fixed orientation. 

Keywords: Parallel kinematic machine. Design, Interval analysis. Comparison, Workspace, Transmission 
factors. 
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1 Introduction 



Parallel kinematic machines (PKM) are known for their high dynamic performances and low positioning errors. 
The kinematic design of PKM has drawn the interest of several researchers. The workspace is usually considered 
as a relevant design criterion [Tl[51[3]. Parallel singularities [3] occur in the workspace where the moving platform 
cannot resist any effort. Thus are very undesirable and generally eliminated by design. The Jacobian matrix, 
which relates the joint rates to the output velocities is generally not constant and not isotropic. Consequently, 
the performances (e.g. maximum speeds, forces, accuracy and stiffness) vary considerably for different points in 
the Cartesian workspace and for different directions at one given point. This is a serious drawback for machining 
applications O [6l [7j . Few parallel mechanisms are isotropic throughout the workspace [H [9j . But their low 
structural stiffness make them inadequate for machining applications because their legs are subject to bending. 

To be of interest for machining applications, a PKM should preserve good workspace properties, that is, 
regular workspace shape and acceptable kinetostatic performances throughout. For example in milling appli- 
cations, the machining conditions must remain constant along the whole tool path [10|, lllj. In many research 
papers, this criterion is not taken into account in the algorithmic methods used to compute the workspace 
volume [121 113j . Other papers present methods that compute the well-conditioned workspace using discretiza- 
tion [14|, I15|. Thus, the results they provide cannot be proved formally. Conversely, interval analysis methods 
applied to well-conditioned workspace computation provide guaranteed results [111 HZ] . 

The comparison of PKM architectures is a difficult but relevant challenge [3 [TH] . Providing tools to allow 
designers or end-users to rigorously compare PKM is indeed necessary since the variety of existing PKM makes 
it hard to choose which one is best-suited for a specific task. 

In this paper, an interval analysis based method is addressed for the design and comparison of 3-DOF 
PKM. This method takes into account two criteria, (i) a regular workspace shape and, (ii) a kinetostatic 
performance index that needs to be as homogeneous as possible throughout the workspace. Two basic tools 
and an algorithm that considers these two criteria are introduced: on the basis of prescribed kinetostatic 
performances, the workspace is analyzed to find out the largest regular dextrous workspace (square, cube, 
cylinder, etc..) enclosed in the Cartesian workspace. 

Two translational parallel mechanisms derived from the Delta robot [2] are compared using this method. 
The first machine, called Orthoglide [19[, features three fixed linear joints which are mounted orthogonally and 
the second one, called UraneSX (Renault Automation) [20j, features three linear joints which arc mounted in 
parallel. In both cases, the mobile platform moves in the Cartesian x — y — z space with fixed orientation. 

Next section presents the interval analysis based method for 3-DOF PKM design. Section 3 presents the 
Orthoglide and UraneSX mechanisms, their kinematic equations and singularity analysis. Section 4 reports 
the comparison between the two mechanisms through the determination of the largest dextrous cube for the 
Orthoglide and the largest dextrous square for the UraneSX enclosed in the workspace. 
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2 Description of the interval analysis based method for 3-DOF trans- 
lational PKM design 



2.1 Preliminaries 

2.1.1 Dextrous Cartesian workspace 

For a 3-axis serial machine-tool, a parallelepiped-shaped Cartesian workspace allows the end-user to visualize 
easily where to place cutting paths. This consideration should also hold for PKM. However the workspace shape 
is often geometrically complex and thus hard to visualize. Therefore, a regular-shaped workspace is needed for 
PKM. Thus, we need to define a regular dextrous workspace which is a regular-shaped workspace included in 
the machine Cartesian workspace. Throughout the dextrous workspace, a kinetostatic performance index (that 
is chosen beforehand) remains as homogeneous as possible. This index can be the local or global conditioning [3j 
of the Jacobian matrix J (that maps the actuated joint rates of the manipulator into the velocity of the mobile 
platform), the force or velocity transmission factors. These last two indices make sense for 3-DOF translational 
PKM with identical actuated joints. 

The method presented in this paper aims at designing such 3-DOF PKM. The velocity transmission factors 
are the ratio between the actuated joints velocities and the velocity of the mobile platform. They are the square 
roots ipi, ip2 and V'a of the real eigenvalues cri, (T2 and of (JJ-^)^^. In order to keep homogeneous kinetostatic 
properties, these factors are bounded inside the dextrous workspace. The regular dextrous Cartesian workspace 
can be defined as a set of points P in the workspace such that ipi, -02 and V'3 are bounded, that is, 

yVoextrous = {P E W \ ^-mm < MP) < Max, « = 1, 2, 3} (1) 

Points P in Woextrous are called dextrous points. 

The values of ipmin and ipmax (resp. Umm and ffmax) depend on given performance requirements. The 
method described further aims at computing the largest dextrous Cartesian workspace included in the Cartesian 
workspace, so that its ratio to the Cartesian workspace is the best one. To be of real interest for milling 
applications, a PKM must indeed include a large regular dextrous workspace in its Cartesian workspace. 

2.1.2 Introduction to ALIAS library 

An algorithm for the definition of the largest dextrous workspace included in the Cartesian workspace is described 
in the following sections. This algorithm uses the ALIAS library [16], which is a C+-I- library of algorithms based 
on interval analysis. These algorithms deal with systems of equations and inequalities whose expressions are an 
arbitrary combination of the most classical mathematical functions (algebraic terms, sine, cosine, log etc..) and 
whose coefficients are real numbers or, in some cases, intervals. An interface exists with Maple that allows the 
automatic generation of C-I--I- codes being given the Maple description of the system and then to compile and 
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run the generated code in order to get the resuh withm the Maple session. Without being exhaustive, ALIAS 
hbrary provides algorithms that enable one to, (i) find an approximation of the real roots of n-dimensional 
systems, (ii) find an approximation of the variety defined by n-dimensional systems, (iii) find an approximation 
of the global minimum or maximum of a function (eventually under equations and/or inequalities constraints) 
up to an accuracy provided by the user, (iv) analyze a system of algebraic equations to determine bounds for 
its real roots. 

2.1.3 Geometric constraints 

The dextrous workspace WDextrous is defined in Eq. [T] as a function of the eigenvalues of (JJ-^)"-'^. These 
eigenvalues are determined by solving the third degree characteristic polynomial V of (JJ-^)^^. To decrease 
the computing time and to avoid numerical problems on singularities, it is recommended to add geometrical 
constraints. These constraints naturally depend on the mechanism architecture (see tj4.2p . 

2.2 A first basic tool: Box verification 

Our purpose is to determine the largest regular dextrous workspace that is enclosed in the Cartesian workspace. 
For a given point, we note valid point if it is a dextrous point and invalid point otherwise. For that purpose we 
need to design first a procedure, called A4{B), that takes as input a Cartesian box B and returns: 

• 1: if every point in B is valid, 

• -1: if no point in B is valid, 

• 0: if neither of the other two conditions could be verified. 

The first step of this procedure consists in considering an arbitrary point of the box (e.g. its center) and to 
compute the eigenvalues at this point: either all of them lie in the range [ffmim '^max] in which case the center 
is called valid or at least one of them lie outside this range and the center point is denoted invalid. 

2.2.1 Valid center point 

In that case if we are able to check that there is no point in B such that one of the eigenvalues at this point 
is equal to (Tmm or amax^ then we can guarantee that every point in B is valid. Indeed assume that at a given 
point B the lowest eigenvalue is lower than cTmin'- this implies that somewhere along the line joining this point 
to the center of the box the lowest eigenvalue is exactly amm ■ 

To perform this check we set the unknown in the characteristic polynomial V of (JJ-^)^^ to amin (and then 
to amax) and we get a polynomial in x, y, z only. We now have to determine if there exists some values for these 
three Cartesian coordinates that cancel the polynomial, being understood that these values have to define a 
point belonging to B. This is done by using an interval analysis algorithm from the ALIAS library [16j. The 
principle of this algorithm is to calculate first the polynomial value for the center point Gb of B. Without 
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lack of generality we may assume that this value is positive. If we are able to determine a point Sb in B such 
that the polynomial value at this point is negative, then we can guarantee that there exists a point on the line 
joining Cb to Sb such that the polynomial is exactly 0. The purpose of the algorithm is now to determine 
if such a point exist. Now let Bi be a box included in B: using interval analysis we are able to calculate 
a range [m^. jAf^.] such that for any point X in Bi we have ruBi < ^{X) < Mb^- Note that this interval 
evaluation is numerically safe as the bounds of the range are calculated by taking into account round-off errors. 
On the other hand these bounds may not be sharp i.e. there may be no X in Bi such that V{X) ~ niBi or 
MBi- Note, however, that the width of the overestimation decreases with the width of Bi. Furthermore we 
may get a sharp evaluation by using, for instance, the derivatives of V. Indeed we may calculate the interval 
evaluation [rx.y.z, Rx.y.z] of dV /dx,y, z and if all three interval evaluations have constant signs (i.e. rj^,y,z > 
or Rx.y.z < 0), then sharp mB^, Mb^ are obtained by setting the variables to fixed values. For instance if > 0, 
then vriBi i-^^Bi) is obtained by setting x to its lower (upper) bound. Note that other methods may also be used 
to determine sharp bounds (see pT ] [22l [23] ) . 
Hence we have the following properties: 

1. if niBi > 0, then for any point in Bi the polynomial V is positive 

2. if MBi < Oj then for any point in Bi the polynomial V is negative 

3. if niBi < and Af^. > and the bounds are sharp, then the polynomial V cancels in Bi 

4. if rriBi < and MBi > and the bounds are not sharp, then we cannot guarantee the sign of the 
polynomial V within Bi 

At that point a simple branch-and-bound algorithm is used: the initial box B is bisected until either all the 
sub-boxes resulting from the bisection satisfy property 1 (in which case we can guarantee that the polynomial 
V never cancels for B and consequently that all the eigenvalues oiV lie in the range [<7min, o'max], which implies 
M{B) = 1) or a sub-box resulting from the bisection satisfies property 2 or 3 which means that at some point 
in B at least one of the eigenvalues of V lies outside the range [amin, o;nax], which corresponds to Ai{B) = 0. 

The algorithm may indeed return for a box that includes only valid points. But the width of this box 
will be lower than a/2 (where a is an accuracy threshold fixed in advance for the computation) and hence the 
final result will be within the tolerance margin of the calculation. The only case in which the calculation will 
be not guaranteed will occur only when a is lower than the machine accuracy. But we may determine that we 
are in such configuration as the width of the box from the machine viewpoint will be 0: if a box of width is 
processed and the algorithm returns 0, then a warning message will be issued indicating that the calculation 
is no more guaranteed. Note, however, that we may still use the algorithm by using a multi-precision package 
such as MPFR that will allow to get a guaranteed result. Furthermore it is doubtful that computing the result 
with an accuracy better than the machine precision makes sense. 
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2.2.2 Invalid center point 

Without lack of generality, we may assume that at the center of the box the largest eigenvalue is greater than 
a„iax ■ If there is no point in B such that one of the eigenvalues is equal to amax , then we can guarantee that for 
any point in B the largest eigenvalue is always greater than (Jmax and consequently A4{B) = — 1. This check is 
performed by using the same method as in the previous case. 

2.3 A second basic tool: Box workspace verification 

During the calculation of the dextrous workspace, we consider a Cartesian box B and we have to examine if this 
box may contain a point that is the center of a Cartesian box Bw with edge length w, which is fully enclosed 
in the robot workspace. We assume here that this workspace is defined by a set of m inequalities Fj such 
that a point X belongs to the workspace if Fj{X) < for all j in [1, m]. Let Bi be a sub-box included in B, 
defined by the three ranges [yi,yi], [zijZi]. All the boxes with edge length w that have as center a point 

in Bi are included in the hull box Hb^ defined by the three ranges [xj — + w/2], [yi — w/2,yi + w/2], 

[zi — 10/2, Zi + iu/2]. As in the previous section we may use interval analysis to compute an interval evaluation 
[mg,,M^,] of all Fj{HBi) with the following properties: 

1. if Mg, < for all j in [1, m], then any point of Bi may be the center of a box with edge length w that is 
included in the workspace 

2. if m^. > 0, then no point of Bi may be the center of a box with edge length w that is included in the 
workspace 

3. if niBi < and M^. > 0, then we cannot determine if some point within Bi may be the center of Bw 

Note also that if the widths of all the ranges defining Bi are lower than w, any box Bw contains the 4 corners 
of the box Bi. 

Using a similar branch-and-bound algorithm as in the previous section, we may now determine if either all, 
none or some points of B may be the center of a box Bw- The initial box B is bisected until either all the 
sub-boxes resulting from the bisection satisfy property 1 (then any point of B may be the center of a box Bw), 
or 2 (no point of B may be the center of a box Bw)- If a sub-box satisfies property 3 and the widths of its 
ranges are lower than w, we check if the corners of B belong to the workspace: if all the corners either belong 
or do not belong to the workspace we continue the bisection. If we have a mixed situation with some corners 
belonging to the workspace whereas other ones do not, we may state that B contains both points that may 
be the center of a box Bw and points that cannot. A similar situation is obtained if we have found at least a 
sub-box that satisfies property 1 and a sub-box that satisfies property 2. 

At that point we may define a procedure G{B,w) that takes as input a box B and an edge length w and 
returns: 

• -1: there is no points in B that may be the center of a box Bw 
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• 1: all the points in B may be the center of a box Bw 

• 0: B contains both points that may be the center of a box Bw and points that cannot. 

2.4 Algorithm for the determination of a cubic dextrous Cartesian workspace 

An algorithm is now described for the determination of a cube that is enclosed in the Cartesian workspace and 
aligned with the coordinate axis, whose edge length is 2w and such that there is no other cube enclosed in the 
workspace with an edge length of 2(w + a). This algorithm can be applied to any 3-DOF manipulator. Other 
shapes for regular dextrous workspace is considered in section [2.51 

The first step is to determine the largest cube enclosed in the workspace with a center located at (0,0,0). 
This is done by using the M procedure on the Cartesian box Binu [~ka, ka], [—ka, ka], [~ka, ka] where k is 
an integer initialized to 1. Each time the A4 procedure returns 1 for Binit (which means that the cube with 
edge length 2ka is enclosed in the dextrous workspace) we double the value of k. If this procedure returns -1 
for a value of k larger than 1 this implies that the cube with edge length ka/2 is in the dextrous workspace 
whereas the cube with edge length ka is not. Hence if fc > 2 (otherwise no improvement is possible) we restart 
the process with k = (fc/2 + fc)/2. After a failure at kfau the principle is to always choose a value of k which 
is the mid-point between the last value k^ of k for which M = 1 and kfaii until kfau = fc^ + 1- For example 
if Ai returns 1 for k=l, 2, 4 and returns -1 for fc = 8 we set fc to 6. Otherwise we have determined that the 
cube with edge length 2fcQ; is enclosed in the dextrous workspace, whereas the cube with edge length 2(fc + l)a 
is not. The value 2fcQ; is hence an initial value for w. Note that the above procedure may be used whatever the 
coordinates of the center: it is implemented as a general purpose procedure C{xmi Um, zm) that takes as input 
the coordinates of a point M and returns the edge length of the largest cube centered at M, that is enclosed in 
the dextrous workspace. 

In the algorithm for determining the largest cube enclosed in the dextrous workspace, we manage a list of 
Cartesian boxes C that are processed by the algorithm in sequence. During the processing, boxes may be added 
to a list. At one step of the algorithm we have n boxes in the list whereas processing box numbered i (which 
means that boxes numbered from to « — 1 have already been processed and may be discarded whereas boxes 
i to n have to be processed). The algorithm stops when all the boxes in C have been processed. The box 
numbered i in the list is denoted Bi and the maximum number of boxes in C is N . 

At the beginning of the algorithm, C has only one box Bq that contains the workspace (for example for the 
Orthoglide Bq= {[— i, L], [— L, L], L, L]}). The algorithm can be described by the following six steps: 

1. calculate w = C(0,0,0) 

2. if i > n EXIT 

3. if G{Bi, w + alpha)= -1, then set z to z + 1 and go to 2 
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4. ii Q [Bi, w + alpha)= 1, then calculate w' = C{xBi,yBi, zb^) where XBi,yBi, ZBi are the coordmates of the 
center of Bi. If w' > w, then update w. Go to step 6 

5. if G{Bi, w + alpha)= 0, then go to step 6 

6. bisect the variable in the box Bi that has the largest range. For example if the box Bi is defined as 
[xijX;'], [yiiVi], [zi,'Zi\ and the variable x has the largest range the bisection process creates two new boxes 
B] = {[xi_, (a^ + x7)/2], [yuvl], [zi_,~i\} and Bj = {[(xi + x7)/2,xl], [yi,W], z*]}- If ?^ < N/2, both boxes 
are stored at the end of the list (and we set otherwise Box B} are stored in C in place of Bi 
whereas box is stored at location i + 1 after a shift of the boxes -Bi+i, . . . , Bn- Set n to rt + 1 and go 
to 2. 

Step 1 allows one to establish an initial value for the maximal edge length. Step 3 eliminates boxes that 
cannot contain the center of the maximal cube due to the workspace limits. Boxes satisfying step 4 are candidate 
to include the center of the largest cube: hence we calculate the largest cube centered at the box that may allow 
to update the current value of the largest edge. Step 6 is the bisection process that allows one to decrease the 
size of the box with the effect of a sharper calculation for the procedure Q. Note also that two storage modes 
that are used for adding the boxes resulting from the bisection process. The second mode allows for a minimal 
memory storage but has the drawback of focusing on a given part of the workspace whereas the center of the 
largest cube may be located in another part. The first mode makes it possible to explore various parts of the 
workspace which may result in large improvement on the value of w but as the drawback of possibly creating a 
large number of boxes. The proposed storage mode allows one to mix the advantages of both storage modes. 

This procedure ensures to determine a cube with edge length w that is enclosed in the workspace and in the 
dextrous workspace, whereas there is no such cube with edge length w -\- a. 

Note that an incremental approach is possible. After having computed w — wi with a given accuracy a it 
is always possible to replace the initial value of w as calculated in step 1 of the algorithm by the value wi when 
computing the cube with a lower value for a. Computation times of the largest cube for various accuracies are 
given for a specific 3-DOF PKM in section [1751 

2.5 Other regular dextrous workspace shapes 

Clearly, considering the largest cube may not be appropriate if the studied PKM has a rectangular or a spherical- 
shaped workspace. The algorithm can thus be modified. Here are for example the necessary changes that must 
be taken into account to consider the largest sphere: the idea is to use spherical coordinates and hence x,y, z 
are substituted by Xc + r sin^p sm9,yc + costl! sin9, Zc + r cos9, with r in [0,R], ip,9 in [0,27r], Xc,yc,Zc being 
the coordinates of the center of the sphere and R its radius. Interval analysis allows to deal with expressions 
involving sine and cosine and hence procedures A4,G can still be used with these new parameters. Similarly 
procedure C{xm ,yM , zm) can be used to determine the largest radius of the sphere centered at {xm, yM, zm) for 
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which the eigenvalues are valid. Hence, with this modification, the algorithm can calculate the largest sphere 
enclosed in the dextrous workspace. 

Spheres and cubes are defined by their center and one additional parameter. But other shapes may involve 
more parameters: for example a cylinder needs a center but also a height and a radius. We can still perform a 
change in the variables so that procedures A4, Q can still be used. The key point is that procedure C(xm, Um, zm) 
has to be modified as we have now two optimization parameters. But in that case volume optimization alone 
has less meaning: for example the optimization result for a cylinder may be a cylinder with a relatively small 
radius and a large height, which may be of no interest. A cylinder of identical radius and height with a lower 
volume may be the most interesting result. A possible way to manage this problem is to assign a range [a, b] 
for the ratio R/h where R is the cylinder radius and h its height. In that case the procedure C has to solve an 
optimization problem which is to maximize the volume of the cylinder under the constraints that the eigenvalues 
are valid and the ratio R/h satisfies a < R/h < b. ALIAS is still able to manage such an optimization procedure. 

2.6 Approximate calculation of the dextrous workspace 

Small modifications in the previous algorithm allow to determine an approximation of the dextrous Cartesian 
workspace Woextrous as a set <S of 3D Cartesian boxes such that for any box B in <S and for any point in 
B the constraints on the eigenvalues are satisfied. The width of all the boxes in the set <S is greater than a 
given threshold e: hence we get only an approximation of the dextrous Cartesian workspace. But the algorithm 
provides the volume 14 of the approximation and a volume error Ve such that the volume Vd of the dextrous 
Cartesian workspace satisfies Vd < Va + Vg. Decreasing the value of e makes it possible to increase Va and to 
decrease Ve- In this paper, this method is used to analyze 3D boxes but it can be applied for any mechanism 
with n d.o.f., the result being a set of nD boxes. 
Initially 14, 14 are set to 0. 

1. if i > n EXIT 

2. if M.{Bi)= -1, then set i to j + 1 and go to 1 

3. if M{Bi)= 1, then store in <S and add its volume to 14 . Set i to i + 1 and go to 1 

4. if M(Bi)= 0, then 

(a) if the largest width of Bi is lower than e, then add its volume to V^, set i to i + 1 and go to 1 

(b) otherwise go to step 5 

5. process bisection for the box B^. Set n to n + 1 and go to 1 

Note that this procedure may be incremental if the boxes neglected at step 4-a are stored in a file !F. Indeed 
a first run with a given e allows to obtain initial values for 14, 14- If the quality of the approximation is not 
satisfactory, we may choose a smaller value of e (e.g. e/2). But instead of starting with the initial Bq, we may 
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use the boxes stored in T , thereby avoiding to repeat computation that has ah'eady been done during the initial 
run. 



3 Description of the Orthoglide and the UraneSX 

The previous interval analysis based design method is now applied to the comparison of two 3-DOF transla- 
tional PKM. It is particularly interesting to compare these two mechanisms because they belong to the same 
architecture family. 

3.1 Orthoglide and UraneSX architectures 

Most existing PKM can be classified into two main families. PKM of the first family have fixed foot points 
and variable length struts and are generally called "hexapods" or "tripods" . PKM of the second family have 
variable foot points and fixed length struts. They are interesting because the actuators are fixed and thus the 
moving masses are lower than in the hexapods and tripods. 

The Orthoglide and the UraneSX mechanisms studied in this paper are 3-DOF translational PKM and 
belong to the second family. Figures [1] and [2] show the general kinematic architecture of the Orthoglide and of 
the UraneSX, respectively. Both mechanisms have three parallel PRPaR identical chains (where P, R and Pa 
stand for Prismatic, Revolute and Parallelogram joint, respectively). The actuated joints are the three linear 
joints. These joints can be actuated by means of linear motors or by conventional rotary motors with ball 




Figure 1: Orthoglide kinematic architecture 



The output body is connected to the linear joints through a set of three parallelograms of equal lengths 
L = AiBi, so that it can move only in translation. Vectors e^ coincide with the direction of the ith linear 
joint. The base points Ai are located at the middle of the first two revolute joints of the i*'* parallelogram, and 
Bi is at the middle of the last two revolute joints of the i*^ parallelogram. 
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For the Orthoglide mechanism, the first hnear joint axis is parallel to the a;-axis, the second one is parallel to 
the y-axis and the third one is parallel to the 2;-axis. When each vector e; is aligned with A^Bi, the Orthoglide 
is in an isotropic configuration and the tool center point P is located at the intersection of the three linear joint 
axes. 




Figure 2: UraneSX kinematic architecture 



The linear joint axes of the UraneSX mechanism are parallel to the z-axis. In fig. [21 points Ai, A2 and A3 
are the vertices of an equilateral triangle whose geometric center is O and such that OAi = R. Thus, points 
i?!, B2 and are the vertices of an equilateral triangle whose geometric center is P, and such that OBi = r. 

3.2 Kinematic equations and singularity analysis 

We recall briefly here the kinematic equations and the singularities of the Orthoglide and of the UraneSX (See 
[201 US] for more details). 

Let 9i and (ii denote the joint angles of the parallelogram about axes ii and ji, respectively (Figs. [Hand 12]) ■ 
Let pi, p2, P3 denote the linear joint variables and L denote the length of the three legs, AiBi. 

For the Orthoglide, the position vector p of the tool center point P is defined in a reference frame (O, x, y, 
z) centered at the intersection of the three linear joint axes (note that the reference frame has been translated 
in Fig. [H for more legibility) . 

For the UraneSX, the position vector p of the tool center point P is defined in a reference frame (O, x, y, 
z) centered at the geometric center of the points Ai, A2, and A^ (same remark as above). 

Let p be referred to as the vector of actuated joint rates and p as the velocity vector of point P: 

P = [Pi h hV, P = [i 2/ zf (2) 
p can be written in three different ways by traversing the three chains AiBiP: 

p = Gipi + {e^ii + X {hi - ai) (3) 
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where and are the position vectors of the points Ai and Bi , respectively, and e.^ is the direction vector of 
the hnear joints, for i ^ 1, 2,3. 

We want to ehminate the three idle joint rates Oi and $i from Eqs. ([3]), which we do by dot-product of 
Eqs. dSl) by b, - a^: 

(b,; - ai)'^p = (b,; - ai)^ei/9i (4) 

Equations Q can now be cast in vector form, namely Ap ~ Bp, where A and B are the paraUel and serial 
Jacobian matrices, respectively: 





(bi 


-air 




Vi 










A = 


(b2 
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V2 





(5) 
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-asr 










V3 





with rji = {hi — aj-^e.; for i = 1, 2, 3. 

Parallel singularities occur when the determinant of the matrix A vanishes, i.e. when det(A) ~ 0. Eq. ^ 
shows that the parallel singularities occur when: 

(bi - ai) = a(b2 - a2) + A(b3 - as) (6) 

that is when the points Ai, A2, B2, A3, and B^, lie in parallel planes. A particular case occurs when the 
links AiBi are parallel: 

(bi--ai) II (b2-a2) and 
(b2-a2) II (bs-aa) and 
(bg-aa) II (bi-ai) 

Serial singularities arise when the serial Jacobian matrix B is no longer invertible i.e. when det(B) = 0. At 
a serial singularity a direction exists along which no Cartesian velocity can be produced. Equation ^ shows 
that det(B) — when for one leg i, {hi — a^) _L e^. 
When B is not singular, we can write, 

p = J^p with J-^ = B-^A (7) 

4 Comparison of the Orthoglide and the UraneSX 

In this section, we calculate the edge length of the largest cube for the Orthoglide, the edge length of the largest 
square for the UraneSX, as well as the location of their respective centers. To simplify the problem, the bounds 
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on the velocity transmission factors are such that tpmin = ^/i^max- 

4.1 Regular dextrous workspace shape 

The Orthoghde and the UraneSX are compared according to the size of their largest regular dextrous Workspace. 
Due to the symmetrical architecture of the Orthoglide, the Cartesian workspace has a fairly regular shape in 
which it is possible to include a cube whose sides are parallel to the planes xy, yz and xz respectively. The 
Cartesian workspace of the UraneSX is the intersection of three cylinders whose axes are parallel to the z- 
axis. Thus, the workspace is unlimited in the z-direction and the Jacobian matrix does not depend on the z 
coordinate. Only the limits on the linear joints define the limits of the Cartesian workspace in the z-direction. 
However, it is possible to include a square in the plane xy. Regular dextrous workspaces are thus chosen to be 
a cube for the Orthoglide and a square for the UraneSX. 

4.2 Geometric constraints 

Section 12.11 is suggested to add geometrical constraints so as to decrease the computing time and to avoid 
numerical problems on singularities. Here, polynomial V is defined only for the points within the intersection 
X of the three cylinders defined by 

x^ +y^ < I? x^ + z'^ < y^ + z^ < (8) 



for the Orthoglide, and, 

{x- R + rf +y'^ < 



a:-iR-r)^y +(y-iR-r)^\ < 



x-{R-r)-j +\^y + {R-r)^j < 

for the UraneSX. With these constraints, matrix B is never singular and thus can be always inverted. To solve 
numerically the above equations and to compare the two mechanisms, the length of the legs is normalized, i.e. 
we set L = 1. 



4.3 Comparison results 

To compare the two mechanisms studied, the leg length L is set to 1 and the bounds on the velocity factor 
amplification are set to ^p ~ [0.5 2], with a = 0.001. For the UraneSX, it is necessary to define two additional 
lengths, r and R. However, the edge length of the workspace depends only on i? — r. 

For the Orthoglide, it is found that the largest cube has its center located at (0.086,0.086,0.086), and that 
the cube edge length is Lworkspace = 0.644. Also, using the incremental approach described in section [2.4) we 
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get for the Orthoglide the computation time of Table [T] on a Sun Blade workstation. 



Accuracy a (mm) 


0.01 


0.001 


0.0001 


0.00001 


Computation time (s) 


360 


150 


504 


900 



Table 1: Computation time of the largest cube enclosed in the dextrous workspace for various 
accuracies 

For the UraneSX, the design parameters are those defined in [20 , which we have normalized to have L — 1, 
i.e. r — 3/26 and R = 7/13. To compare the two mechanisms, we increase the value of R such that R' = R + X 
with A — [0.0, 0.2]. For R < 7/13, the constraints on the velocity amplification factors are not satisfied. 



A 


Center 


-^Workspace 


0.00 


(-0.0178,-0.0045) 


0.510 


0.05 


(-0.0179,-0.0022) 


0.470 


0.10 


(-0.0225,-0.0031) 


0.420 


0.15 


(-0.0245,-0.0018) 


0.370 


0.20 


(-0.0211,-0.0033) 


0.320 



Table 2: Variations of the edge length of the square workspace for the UraneSX mechanism 

The optimal value of R' is obtained for A = 0, i.e. for the design parameters defined in [20 for an industrial 
application (see table [2]). To expand this square workspace in the z-direction, the range limits must be equal 
to the edge length of the square plus the range variations necessary to move throughout the square in the x ~ y 
plane. 

The constraints on the velocity amplification factors used for the design of the Orthoglide are close to those 
used for the design of the UraneSX which is an industrial machine tool. For the same length of the legs, the 
size of the cubic workspace is larger for the Orthoglide than for the UraneSX. 

For the Orthoglide, the optimization puts the serial and parallel singularities far away from the Cartesian 
workspace [T^ . The UraneSX has no parallel singularities due to the design parameters [R — r < L), but 
serial singularities cannot be avoided with the previous optimization function. To produce the motion in the 
z-direction, the range limits of the linear joints are set such that the constraints on the velocity amplification 
factors are not satisfied throughout the Cartesian workspace. 

The range limits Api of each prismatic joint can be decomposed into two parts. For the Orthoglide (resp. 
for the UraneSX), the first part Afi makes it possible to move the mobile platform throughout the face of the 
prescribed cube that is perpendicular to the considered prismatic joint axis (resp. throughout the prescribed 
square). The second part is equal to the edge length of the cubic workspace Lworkspace- The equations of the 
inverse kinematic model allow us to compute Afi for the two mechanisms. 

For the Orthoglide, the position and the size of the prescribed cube define three range limits for the x — y — z 
platform coordinates, 

X = [-0.322 + 0.085,0.322 + 0.085] (9a) 
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y = [-0.322 + 0.085,0.322 + 0.085] (9b) 
z = [-0.322 + 0.085,0.322 + 0.085], (9c) 

For the UraneSX, the position and the size of the prescribed square define two range limits for the x — y platform 
coordinates, 

X =[-0.255- 0.018,0.255- 0.018] (9d) 
y =[-0.255,0.255]. (9e) 

For the Orthoglide, all A/; are equal due to the symmetrical architecture. For the UraneSX, we take A/ = 
Max{Afi). The results are A/ = 0.181 and Ap = 0.825 for the Orthoglide and A/ = 0.353 and Ap = 0.863 
for the UraneSX. This means that the range limits are quite similar for the same leg length. To calculate the 
volume of the Cartesian workspace of the two mechanisms for the previous range limits, we have used a CAD 
system. Results are given in table |31 





Cartesian workspace volume 
with optimized ranges limits 


cubic dextrous workspace volume 


ratio 


Orthoglide 


0.566 


0.265 


46.8% 


UraneSX 


0.544 


0.132 


24.3% 



Table 3: Workspace volumes of the two mechanisms 



To help understand these results. Fig. [3] and H] show the location of the largest cubic workspace inside the 
Cartesian workspace. As the Cartesian workspace of the Orthoglide is regular and admits a quasi-cubic shape, 
the ratio between the cubic workspace and the Cartesian workspace is better than for the UraneSX. 




Figure 3: Cartesian workspace and dextrous workspace for the Orthoglide mechanism with opti- 
mized range limits 

In tableUl the design parameters are compared to achieve the same cubic dextrous workspace with Lworkspace = 1- 
The legs length is directly connected to the dynamic properties of the mechanism. The range limits and the 
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Figure 4: Cartesian workspace and dextrous workspace for the UraneSX mechanism with opti- 
mized range limits 

legs length are important parameters in the determination of the total size of the mechanism and in its global 
cost. The volume of the Cartesian workspace allows us to characterize the shape and the volume of motion of 
the tool with regard to the useful Cartesian workspace dedicated to manufacturing tasks (cubic workspace). 





Leg length 


Range limits 


Volume of the Cartesian workspace 


Orthoglide 


1.55 


1.28 


2.13 


UraneSX 


1.96 


1.69 


4.12 



Table 4: Synthesis of the comparative study for the same cubic Cartesian workspace 

These criteria allow us to optimize some geometric parameters to design a machine tool for milling appli- 
cations. Although in this approach, the kinetostatic properties of the Orthoglide are better than the UraneSX 
ones, we cannot assert that the Orthoglide is better than the UraneSX. One reason is that these two PKM 
are not aimed at identical manufacturing tasks. The main applications of the UraneSX are drilling, facing and 
tapping whereas the Orthoglide is more universal. 

Other shapes of regular dextrous workspaces can be computed for the Orthoglide and the UraneSX by using 
cylindrical or spherical coordinates to have the largest cylinder or sphere respectively, even if these shapes are 
generally less relevant for milling applications. 

5 Conclusions 

This paper introduces an interval analysis based study for the design and the comparison of 3-DOF PKM. Two 
basic tools and an algorithm are described to determine the largest regular dextrous workspace enclosed in the 
Cartesian workspace. The dextrous workspace is a part of the Cartesian workspace in which the velocity ampli- 
fication factors remain within a predefined range. This means that throughout the dextrous workspace, milling 
tool paths arc available because the variations of the kinematic performances index remain under reasonable 
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values. The regular dextrous workspace shape is a cube for the Orthoglide and a square for the UraneSX. 
This general method is coupled with geometric constraints associated with the mechanisms studied to avoid 
numerical problems at singular configurations. The shape of the dextrous workspace was chosen for milling 
applications but it can be different for other applications. The range limits and the volume of the Cartesian 

workspace were calculated to compare the two mechanisms. 
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